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PREFACE 


The  work  reported  herein  was  conducted  at  the  Arnold  Engineering 
Development  Center  (AEDC),  Air  Force  Systems  Command  (AFSC), 
under  the  sponsorship  of  Air  Force  Cambridge  Research  Laboratories 
(AFCRL),  AFSC,  under  Program  Element  65802F.  The  AFCRL  proj¬ 
ect  monitor  was  R.  E.  Good,  who  provided  the  technical  objectives  of 
the  project.  The  results  were  obtained  by  ARO,  Inc.  {a  subsidiary  of 
Sverdrup  &  Parcel  and  Associates,  Inc.),  contract  operator  of  AEDC, 
AFSC,  Arnold  Air  Force  Station,  Tennessee.  This  work  was  conducted 
from  November  1972  to  July  1973  under  ARO  Project  No.  VF233,  and 
the  manuscript  (ARO  Control  No.  ARO- VKF-TR- 74- 13)  was  submitted 
for  publication  on  January  18,  1974. 


1 


A  E  D  C-T  R  -74-32 


CONTENTS 


Page 


1.0  INTRODUCTION  AND  THEORY . 5 

2.0  EXPERIMENTAL  APPARATUS  AND  METHOD . 14 

3.0  EXPERIMENTAL  RESULTS  AND  ANALYSIS . 16 

4.  0  DISCUSSION  AND  CONCLUSIONS . 32 

REFERENCES  .  .....  34 


ILLUSTRATIONS 


Figure 

1.  Experimental  Configuration . '8 

2.  Filter  Transmission  Curve . 15 

3.  Condensation  Diagram  of  Argon  for  a  Sonic 

Orifice  Expansion . 17 

4.  Rayleigh  Scattering  Results  for  PQ  =  200  torr . 18 

5.  Rayleigh  Scattering  Results  for  PQ  =  250  torr . 18 

6.  Rayleigh  Scattering  Results  far  PQ  =  350  torr . 19 

7.  Rayleigh  Scattering  Results  for  PQ  =  500  torr . 19 

8.  Rayleigh  Scattering  Results  for  PG  =  750  torr . 20 

9.  Axial  Variation  of  Effective  Argon  Density  (N/N0)eff  ...  20 

10.  Pressure  Variation  of  Condensation  Onset  Location 

of  Argon . 22 

11.  Pressure  Variation  of  Super  saturation  Ratio . 23 

12.  Pressure  Variation  of  Degrees  of  Supercooling . 23 

13.  Axial  Variation  of  Depolarization  Ratio  of  Argon . 24 

14.  Axial  Variation  of  Ratio  R . 25 

15.  Variation  of  (N/Z)*  with  PQ . 27 

16.  Scaled  Variation  of  (N/Z)*  with  PQ . 28 


3 


AE  DC-TR-74-32 


Figure  Page 

17.  Poisson  Distribution  Function  for  a  Mean  Value  of  20  .  .  .  29 

18.  Variation  of  Mass  Fraction  g  with  Reservoir 

Pressure  PQ . 30 

TABLES 

1.  Condensate  Mass  Fraction  Results . 25 

2.  Cluster  Size  Results  Using  Mass  Fraction  g  of  Ref.  5  ....  29 

NOMENCLATURE . 36 


4 


A  E  D  C-TR-74-32 


1.0  INTRODUCTION  AND  THEORY 


1.1  INTRODUCTION 

Interest  in  the  understanding  of  the  fluid  mechanics  and  rate  kinet¬ 
ics  of  two-phase  flow  fields  results  not  only  from  a  desire  to  exploit 
such  flows  for  the  study  of  dimer  and  cluster  collision  processes,  but 
also  from  a  preference  for  avoiding  their  occurrence  by  judicious  se¬ 
lection  of  the  reservoir  parameters,  if  possible.  In  the  latter  category 
are  included  aerospace  simulation  facilities,  where  such  processes  af¬ 
fect  the  aerodynamic  simulation  parameters,  and  fluid  mechanical  de¬ 
vices,  such  as  turbines,  'which  suffer  loss  in  efficiency  as  a  result  of 
condensation.  The  authors’  interest  in  condensing  supersonic  and  hyper¬ 
sonic  flow  fields  arises  from  a  desire  to  characterize  further  the  cluster 
growth  process  resulting  from  homogeneous  condensation  for  the  purpose 
of  understanding  both  the  dissimilar  areas  of  chemical  kinetics  in  two- 
phase  flows  and  enhanced  scattering  by  the  clusters  of  radiation  inci¬ 
dent  upon  a  condensed  flow  field. 

Several  previous  studies  of  condensation  in  expanding  gas  flows 
have  considered  the  cluster  growth  of  argon  (Ar)  because  it  is  a  non¬ 
toxic,  nonreactive,  inexpensive,  noble  gas  atom  which  is  believed  to 
be  capable  of  description  using  a  spherically  symmetric  Lennard- Jones 
12-6  interatomic  potential.  Such  spherical  symmetry  greatly  simplifies 
theoretical  calculations  of  equilibrium  Ar2  dimer  concentrations,  the 
values  of  which  are  required  for  rate  calculations.  Additionally,  the 
absence  of  internal,  ground  electronic  state  modes  of  motion  precludes 
the  existence  of  relaxation  phenomena  in  the  flow  field  and  enables  the 
accurate  prediction  of  gas  dynamic  properties  throughout  the  flow  field. 
Finally,  and  of  prime  importance,  Ar  condenses  readily.  The  majority 
of  these  previous  studies  (Refs.  1  and  2)  have  relied  on  mass  spectro- 
metric  detection  of  the  clusters,  or  polymers,  and  despite  the  inherent 
difficulties  of  the  uncertainty  of  ionization  cross  sections  and  diffusive 
separation  effects,  this  method  of  measurement  has  definitively  indi¬ 
cated  the  onset  of  condensation.  It  should  be  noted  that  such  detection 
methods,  although  not  all  alike  in  sampling  procedures,  have  the  com¬ 
mon  feature  of  being  located  far  from  both  the  flow  source  and  the  point 
of  condensation  onset.  As  a  result,  the  detected  polymer  signals  are 
the  result  of  integration  of  the  rate  equations  along  the  streamtube  of  the 
flow  through  normally  large  gradients  of  both  temperature  and  density, 
a  most  undesirable  effect.  Consequently,  one  must  resort  to  assuming 
various  kinetic  models,  performing  the  required  integrations  and  com¬ 
paring  the  result  with  the  experimental  values  (Ref.  3). 
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Recently,  electron  diffraction  studies  (Ref.  4)  have  been  performed 
on  argon  cluster  beams,  and,  although  these  measurements  have  indi¬ 
cated  the  existence  of  large  clusters,  quantitative  estimates  of  the  con¬ 
centration  of  the  clusters  have  not  been  possible.  It  should  be  noted 
that  the  experimental  configuration  of  observing  a  sampled,  or  skimmed, 
beam  is  used  for  such  diffraction  studies  just  as  it  is  for  mass  spectro- 
metric  detection;  consequently,  considerations  concerning  the  observa¬ 
tion  of  an  integrated  effect  apply  equally  to  the  diffraction  method  of 
observation  so  long  as  collisionless  sampling  methods  are  employed. 

The  work  of  Audit  (Ref.  4)  has  demonstrated  the  utility  of  the  technique 
for  distinguishing  between  liquid  and  solid  phase  condensation  by  obser¬ 
vation  of  the  diffraction  pattern  and  the  spatial  extent  of  ordering  present. 
Moreover,  estimates  of  interatomic  spacing  and  temperature  of  the  crys¬ 
talline  structure  are  both  possible,  and  this  information  is  most  useful 
for  determining  the  partition  of  enthalpy  between  the  condensate  and  gas 
phases  of  the  flow  field.  Recently,  measurements  (Ref.  5)  of  the  mono¬ 
mer  flow  speed  using  mass  spectrometric  methods  have  enabled  the  fi¬ 
nal  mass  fraction  of  condensate  to  be  determined,  and  the  electron  beam 
fluorescence  method  (Ref.  6)  has  been  applied  to  the  task  of  determining 
condensate  mass  fractions  for  the  far  field,  low-density  region  of  nitro¬ 
gen  (N2)  condensing  flows.  Moreover,  Rayleigh  scattering  measure¬ 
ments  (Ref.  6)  are  capable  of  determining  within  the  flow  field  the  spatial 
location  of  condensation  onset  as  well  as  following  the  cluster  growth 
process.  Additionally,  in  the  early  stages  of  cluster  growth  when  only 
dimers  exist,  quantitative  estimates  of  dimer  mole  fraction  are  pro¬ 
vided  by  such  measurements.  This  information,  of  course,  provides 
a  rather  stringent  test  for  condensation  theories,  and  the  test  itself 
does  not  suffer  such  possible  ambiguities  as  sampling  perturbations 
and  cluster  dissociation  by  electron  impact  as  does  the  test  of  the  mass 
spectrometer. 

Historically,  the  observation  of  dimers  of  Ar  in  a  sampled,  molec¬ 
ular  flow  by  Leckenby  and  Robbins  (Ref.  7)  signaled  the  beginning  of  re¬ 
cent  intensive  investigations  of  expansion  flow  fields  by  Golomb,  et  al. 
(Refs.  1  and  2),  Hagena  and  coworkers  (Ref.  8  and  9),  and  Milne  et  al. 
(Ref.  3).  It  should  be  noted  that  the  authors  in  Refs.  1  through  3  all 
used  conventional  quadrupole  or  magnetic  sector  mass  spectrometers 
for  the  far-field  detection,  whereas  Hagena  employed  a  specially  con¬ 
structed  retarding  field  device  capable  of  determining  a  quantity  resem¬ 
bling  the  mean  value  of  N/Z,  the  ratio  of  number  of  atoms  per  cluster 
and  the  charge  of  the  cluster. 
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Since  the  primary  interest  at  present  is  in  the  nonequilibrium  pro¬ 
duction  of  condensate  in  expansion  flows,  the  notable  work  of  Leckenby 
and  Robbins  is  of  importance  to  this  work  only  for  historical  purposes, 
and  attention  to  previous  work  will  be  directed  to  that  of  Refs.  1  and  2 
as  well  as  the  most  recent  work  of  Yealland  et  al.  (Ref.  10).  It  is  in¬ 
teresting  to  note  the  varying  degrees  of  confidence  by  these  investiga¬ 
tors  in  quantitative  interpretation  of  their  respective  data.  The  work  of 
Golomb  et  al.  (Ref.  2)  is  noteworthy  for  attention  to  the  problems  and 
ambiguities  arising  from  the  sampling  process  itself,  which  occur  in 
density  regions  commonly  considered  as  free  molecular.  Definite  dif¬ 
ferences  in  results  are  noted,  for  example,  whether  the  skimmer  is 
cooled  or  uncooled  and  whether  the  mass  spectrometer  extracts  the 
ions  either  parallel  or  perpendicular  to  the  flow  velocity.  In  light  of 
such  differences,  Golomb  et  al.  do  not  place  great  reliance  upon  the 
absolute  values  of  their  results  but  do  contend  that  the  onset  of  massive 
condensation  is  observed.  In  addition  to  observation  of  the  increasing 
concentration  of  polymers  in  the  sampled  flow,  the  study  of  Golomb  et 
al.  also  considered  variations  in  measured  flow  speed  with  pressure  as 
an  indication  of  condensation  onset,  and  Bailey  (Ref.  5)  has  determined 
the  mass  fraction  of  condensate  of  various  species  using  this  technique. 

An  opposing  view  of  the  confidence  to  be  given  quantitative  inter¬ 
pretation  of  mass  spectrometric  data  is  obviously  held  by  Milne  et  al. 
who  go  so  far  as  to  use  their  results  to  attempt  to  determine  the  func¬ 
tional  dependence  of  the  various  rate  constants.  With  similar,  but  less¬ 
er,  fervency,  Yealland  et  al.  present  quantitative  functions  for  the 
prediction  of  dimer  mole  fractions  for  the  expansion  flow  fields  they  in¬ 
vestigated.  Hagena  asserts  his  results  are  quantitatively  accurate  to 
the  extent  that  he  can  characterize  the  expansion  by  some  type  of  mean 
value  of  cluster  size.  Finally,  authors  in  both  Refs.  2  and  8  have  con¬ 
sidered  in  detail  the  scaling  of  condensation  onset  or  cluster  size  with 
orifice  diameter,  reservoir  pressure,  temperature,  and  gas  specie. 

A  serious  difficulty  encountered  in  comparing  results  of  these  in¬ 
vestigations  is  the  necessity  of  extrapolating  various  parameters  based 
upon  the  experimentally  determined  scaling  laws.  Although  such  an  ex¬ 
trapolation  is  questionable  and  undesirable,  no  other  choice  exists.  As 
will  be  discussed  in  detail  in  a  later  section,  such  a  comparison  of  the 
results  of  Refs.  2  and  8  does  not  yield  agreement.  For  example,  for 
reservoir  pressures,  PQ,  for  which  Golomb  et  al.  observed  a  prepon¬ 
derance  of  dimers  of  Ar2,  Hagena  and  Obert  predict  the  mean  cluster 
size  to  be  20.  The  mean  cluster  size  of  Ar,  (N/Z)*,  scales  with  ori¬ 
fice  diameter,  D,  as  P0D®*  ^  according  to  Ref.  8,  whereas  the  dimer 
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peak  of  Ref.  2  scales  essentially  as  Pod0-  5 1  and  the  situation  is  no 
better  with  regard  to  scaling  laws  relating  PQ  and  reservoir  tempera¬ 
ture,  T0.  The  purpose  of  this  work  was  to  obtain  information  concern¬ 
ing  the  spatial  location  of  condensation  onset,  supersaturation  ratios 
as  a  function  of  reservoir  pressure,  and  cluster  growth  processes,  in¬ 
cluding  estimates,  if  possible,  of  the  size  of  the  scattering  cluster. 


1.2  RAYLEIGH  SCATTERING 

Since  thorough  discussions  of  photon-atom  and  molecule  scattering 
processes  have  been  presented  elsewhere  (Ref.  11),  only  a  brief  review 
of  the  phenomenon  will  be  considered  here.  Attention  will  be  restricted 
to  photon-molecule  scattering  processes  for  which  no  change  in  internal 
energy  states  exists,  and  this  process  will  be  designated  as  Rayleigh; 
scattering.  The  laser  source  of  intensity,  IQ,  and  wavelength,  X,  pro¬ 
vides  a  photon  beam  perpendicular  to  both  the  gas  flow  velocity  and  the 
unit  vector  in  the  direction  of  the  collection  optics  for  the  scattered  ra¬ 
diation.  Figure  1  shows  this  configuration  as  well  as  the  polarization 


+z,  Scattering  Direction 


Figure  1.  Experimental  configuration. 
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vectors  of  the  incident  and  scattered  radiation.  It  is  known  that  inten¬ 
sity,  Ig,  with  wavelength,  A,  scattered  from  a  collection  of  uniform 
scattering  sources  of  number  density,  N,  is  given  by 


'A  -  KaNa2/A* 


(i) 


where  Ka  includes  various  transmission  and  optical  sensitivity  constants 
as  well  as  fundamental  atomic  constants  and  a  is  the  polarizability  of  the 
scatterer. 


The  polarizability  of  the  monomer  Ar  is  known  (Ref.  11)  to  be  1.64 
A^.  Since  spinless,  ground- state  Ar  atoms  of  spherical  symmetry 
will  be  considered,  the  depolarization  ratio  is  zero  for  gas  densities 
sufficiently  low  that  insignificant  numbers  of  quasi- molecules  Ar2 
exist. 


Following  onset  of  condensation,  a  collection  of  clusters  or  poly¬ 
mers  of  various  sizes  exist  within  the  flow,  and  the  number  density  of 
the  j-mer  of  the  collection,  i.e.,  Arj,  is  denoted  by  Nj;  the  polariza¬ 
bility  of  the  j-mer  is  designated  as  nj.  The  scattering  equation  for 
such  an  assembly  is  given  by 


‘A-  >  x- 


OO  0 

2  N.a? 


i=l 


(2) 


or 


where 


>A  -  ayAVi 


+  2  (N./NT)(a./a1)' 


T/  i=2  1 


NT  =  2  N. 
T  i=l  1 


(3) 


Using  the  reservoir  number  density  for  nondimensionalizing  the  im¬ 
portant  parts  of  the  equations, 


IA  =  x  (Nx/No)  x 


+  iw  x  fli  1 


(4) 
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For  isentropic,  noncondensing  expansions  of  the  monomer  from 
the  same  reservoir  density,  N0,  Eq.  (4)  Would  be 

<w°  -  (5) 

where  the  form  of  (N^/Nq)0  is  known  (Ref,  12)  as  a  function  of  axial 
distance,  x/D,  and  specific  heat  ratio,  y. 

The  ratio  of  Eq.  (4)  and  (5)  yields 

[(Is/Io)/(Is/Io)0]  -  (Nj/Nj)  =  (NT/Np  x  Xc  x  S  P(i)(a./a])2  (6) 

i=  2 

where  the  mole  fraction  of  condensate,  Xc,  is  defined  a$ 

oo 

Xc  =  2  Ni/NT  =  Nc/NT  <7) 

and  the  probability  of  the  existence  of  an  i-mer  p(i)  is  defined  as 

/oo 

IN,  (8) 

The  condensate  mass  fraction,  g,  is  defined  as 

g  =  I  KN./N T)/f  i(N./NT)  (9) 

i=2  1  1  /  i=l  11 

or 

OO  CO 

g  =  (Nc/Nt)  X  Ip(i)  X  i/E(N1/NT>  +  (Nc/NT)S2ip(i)]  (10) 

For  the  purpose  of  data  interpretation  and  comparisons,  Eqs.  (6)  and 
(10)  will  be  specialized  to  the  following  three  cases: 

1.  The  cluster  .distribution  is  characterized  by  a  Dirac  delta 
function  such  that  all  clusters  are  dimers,  and  the  mole 
fraction  of  condensate  is  much  less  than  unity, 

2.  The  cluster  distribution  is  characterized  by  a  Dirac  delta 
function  such  that  all  clusters  are  J-mers,  where  J  »  1, 
and  the  mole  fraction  of  condensate  is  much  less  than  unity,, 
and 

3.  The  cluster  distribution  is  characterized  by  a  Poisson  dis¬ 
tribution  function  with  mean  cluster  size  of  J  »  1. 


10 


AEDC-TR-74-32 


For  the  first  case,  since  the  mole  fraction  of  condensate  Xc  «  1, 

N°  a  Nj  «  Nt 


so  Eqs.  (6)  and  (10)  become 


R  =  ^2  x  (a2/aj)2  (11) 

and 

g  =  2X2/(1  +  X2)  (12) 

where  X^  is  the  mole  fraction  of  the  i-mer  and 

R  =  (I  /I  )/(I  /I  )°  -  1 

S  O  S  0 

Since  the  Ar2  dimer  is  a  weakly  bound,  van  der  Waals-type  molecule, 
it  is  reasonable  to  assume  (Ref.  11) 

CL  2  ^  CL  j  as  2 


Therefore, 


R  -  4X2  (13) 

Consequently,  for  expansion  flow  fields  for  which  g  «  1  and  only  dimers 
exist  throughout  the  flow  field,  a  measurement  of  the  scattering  function 
R  provides  a  direct  measure  of  X2  by  means  of  Eq.  (13),  and,  of  course, 
g  is  determined  using  X2  and  Eq.  (12).  Although  the  reservoir  condi¬ 
tions  may  be  such  as  to  produce  a  flow  field  containing  sufficient  con¬ 
densate  in  the  far  field  to  invalidate  the  assumptions  of  Case  1,  a  region 
remains  in  the  spatial  vicinity  of  condensation  onset  for  which  the  as¬ 
sumptions  of  Case  1  are  valid.  Intuitively,  one  expects  such  a  region 
of  validity  to  become  more  narrow  in  spatial  extent  as  the  reservoir 
conditions  are  altered  in  the  direction  of  enhancement  of  condensation. 

Considering  Case  2  for  which  p(i)  equals  unity  for  i  =  J  and  zero 
otherwise  and  Xc  «  1,  Eq.  (6)  becomes 

R  «  Xc(aj/a,)2  (14) 

where  it  is  assumed  that 


= 


Nj  =  N,r 
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If  it  is  assumed  that 


aj/ai  =  J 


Eq.  (14)  becomes 


R  =  J2Xc  (15) 

It  should  be  noted  this  approximation  regarding  the  polarizability  of  the 
j-mer,  J  »  1,  becomes  invalid  when  long-range  order  or  collective  ef¬ 
fects  are  significant.  The  upper  limit  in  J  for  which  this  approximation 
is  valid  is  not  known  at  present.  The  mass  fraction  of  condensate  is 
given  by 


g  =  JXc/[l  +  (J  -  l)Xc]  (16) 

Considering  Case  3  for  which  the  distribution  function  is  character¬ 
ized  by  a  Poisson  distribution  of  mean  cluster  size  of  J 

P(i)  =  e-J  JVi!  (17) 

It  is  assumed  that  only  the  clusters  are  described  by  such  a  distribution; 
the  monomer  is  not  included.  As  will  be  shown  later,  the  distributions 
to  be  considered  will  be  such  that  the  mean  cluster  size  J  will  greatly 
exceed  unity  and  the  contribution  to  the  usual  Poisson  distribution  by 
p(l)  is  negligible.  That  is  to  say,  because  of  the  large  values  of  J,  the 
normalization  of  the  Poisson  distribution  function  is  not  significantly 
affected  by  our  neglect  of  p(l). 

As  a  first  approximation,  it  is  assumed  that  Ni/N'J  and  are 

approximately  unity,  so  Eq.  (6)  is  written  as 

R  =  X  2  p(i)(a./a1)2  (18) 

ci— 2  1 

where  p(i)  is  given  by  Eq.  (17).  Concerning  the  approximation 

Nj/N°  =  1 

if  such  is  not  the  case  it  is  probable,  and  experimentally  observed,  that 

w<w> »  i 
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so  that  little  error  is  encountered  on  the  left-hand  side  of  Eq.  (6)  by 
assuming  Nj  «  Nj.  The  primary  error  is  expected  to  arise  on  the  right- 
hand  side  of  Eq.  (6)  in  approximating  N^p  by  Nj.  However,  it  is  sug¬ 
gested  that  any  error  introduced  by  such  an  approximation  is  perhaps 
insignificant  for  the  flow  fields  investigated  in  this  study. 

Once  again  approximating  a-J a\  by  i,  it  is  seen  that  the  second  mo¬ 
ment  of  the  Poisson  distribution  is  required  for  the  evaluation  of  R  in 
Eq.  (18).  It  is  easily  shown  that 

I  i2P(i)  =  J2  +  J 

i-2 

where 

oo 

J  =  2  ip(i) 

i=2 

and  is  the  mean  cluster  size  of  the  distribution.  Therefore, 

R  -  Xc(J2+J)  (19) 

and  g  is  given  by 

g  =  XCJ/[1  +  XC(J  -  1)]  (2  0) 

or  for  J  »  1 

g  =  XCJ/(1  +  XCJ).  (21) 

The  approach  taken  for  the  analysis  of  the  experimental  results  will 
be  to  obtain  experimental  values  of  R,  use  the  results  of  Ref.  8  for  val¬ 
ues  of  J,  calculate  Xc  and  g,  and  compare  these  results  with  those  of 
Golomb  et  al.  (Ref.  2)  and  Bailey  (Ref.  5). 

Finally,  it  should  be  noted  that  measurement  of  the  depolarization 
ratio,  p,  through  the  cluster  growth  region  provides  an  interesting  ac¬ 
count  of  the  change  of  spatial  symmetry  of  the  scatterer.  In  the  event 
of  scattering  by  only  monomers,  p  of  Ar  is  zero.  As  the  cluster  growth 
proceeds  through  the  dimer  and  trimer  stages,  one  would  intuitively  ex¬ 
pect  a  reduction  in  degree  of  symmetry  and  an  increase  from  zero  of 
the  depolarization  ratio.  As  the  cluster  grows  and  assumes  larger  di¬ 
mensions,  it  is  expected  that  the  degree  of  symmetry  should  once  again 
increase  and  p  should  decrease.  Consequently,  a  qualitative  picture  of 
the  growth  process  will  result. 


13 


AEDC-TR-74-32 


2.0  EXPERIMENTAL  APPARATUS  AND  METHOD 


2.1  GAS  SOURCE  AND  VACUUM  CHAMBER 

Since  the  experimental  apparatus  is  essentially  the  same  as  that 
described  in  Ref.  6,  only  a  brief  summary  of  the  system  will  be  pre¬ 
sented. 

The  sonic  orifice  source  was  of  3.  2-mm  diameter  with  a  diameter- 
to-thickness  ratio  of  approximately  25.  The  orifice  was  located  at  the 
end  of  a  tubular  stainless  steel  reservoir  of  1-cm  inner  diameter,  and 
equilibration  of  the  gas  and  tube  wall  temperatures  was  assured.  Res¬ 
ervoir  temperature  and  pressure  were  measured  using  calibrated  trans¬ 
ducers,  and  the  flow-field  data  were  acquired  using  argon  of  both  high 
and  ultrahigh  purity,  the  latter  being  99.  999  percent  pure,  and  no  differ¬ 
ence  in  the  results  was  noted.  Two  0.  025-/um  particulate  filters  were 
installed  in  the  inlet  gas  line  to  minimize  heterogeneous  condensation 
processes.  The  gas  source  was  mounted  on  a  traversing  mechanism 
to  provide  three  degrees  of  freedom  and  enable  flow-field  profile  meas¬ 
urements  with  fixed  optical  instrumentation.  The  movement  mechanism 
has  an  accuracy  and  reproducibility  of  0.  013  cm  in  the  axial  direction, 
and  a  more  detailed  description  and  sketches  are  given  in  Ref.  6. 

The  vacuum  chamber  enclosing  the  motor- driven  traversing  mecha¬ 
nism  was  cryopumped  using  liquid  N2  and  gaseous  helium  (He),  and 
background  pressures  of  lO-^  and  less  than  10“^  torr  were  attained  for 
no  gas  flow  and  maximum  mass  flow  rate,  respectively.  Pressure  and 
temperature  fluctuations  in  the  reservoir  were  on  the  order  of  one  to 
two  percent  in  the  worst  case  and  typically  were  much  smaller.  A  more 
complete  description  of  both  the  gas  dynamic  instrumentation  and  cham¬ 
ber  systems  is  found  in  Ref.  6. 

2.2  RAYLEIGH  SCATTERING  APPARATUS 

Figure  1  shows  the  experimental  configuration  employed  for  the 
laser  scattering  studies.  A  Coherent  Radiation  Laboratories  Modelo 
52-B  argon  ion  laser  was  used  with  incident  wavelength  set  at  4880  A 
and  a  0.  5-watt  intensity.  A  half-wave  plate  was  used  on  the  input  beam 
to  rotate  the  polarization  vector  by  90  deg,  and  HN-22  Polaroid®  mate¬ 
rial  was  inserted  in  the  detection  train  to  select  either  the  parallel  or 
perpendicular  polarization  component  of  the  scattered  radiation.  The 
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rejection  ratio  for  transmission  of  crossed  polarization  vector  intensity 
was  105.  Details  of  the  optical  apertures  used  to  reduce  background 
scattering  are  presented  in  Ref.  6  and  will  not  be  repeated  here.  An 
International  Telephone  and  Telegraph  (ITT)  FW130  photomultiplier 
tube  cooled  to  -26°C  was  used  for  the  photon  counting  detection  method. 
For  pulse  processing  and  registration,  an  Ortec  system  consisting  of 
amplifier,  discriminator,  and  either  a  dual  counter/ timer  or  a  rate- 
meter  was  utilized.  Spectral  isolation  was  accomplished  by  means  of 
a  10-A  half-width  transmission  filter  centered  at  4880  A,  the  trans¬ 
mission  curve  of  which  is  shown  in  Fig.  2.  The  sample  volume  was 
on  the  order  of  1  mm^.  Linearity  of  the  output  signal  with  incident 
laser  power  was  verified.  The  pressure  variation  of  the  Rayleigh 
scattered  radiation  from  Ar  has  been  measured  over  three  orders  of 
magnitude  range  in  pressure,  and  linearity  has  been  observed.  Fi¬ 
nally,  room  temperature  measurements  of  the  depolarization  ratio 
of  Ar  were  performed  for  several  pressures  and  it  was  found  that 
p^lx  10~4. 


Figure  2.  Filter  transmission  curve. 


15 


AEDC-TR-74-32 


The  procedure  for  alignment  of  the  laser  beam  with  respect  to  the 
gas  source  in  the  axial  direction  was  accomplished  using  a  pitot  probe 
of  fixed  location  within  the  chamber,  and  Ref.  6  describes  this  procedure. 
The  centerline  of  the  flow  was  found  using  radial  profiles  and  selecting 
the  center  of  symmetry.  Such  profile  measurements  were  required  for 
several  axial  positions  to  define  the  direction  of  the  axis  of  the  cylindri- 
cally  symmetrical  flow  field. 

Although  calibrations  of  the  constant  Ka  in  Eq.  (1)  were  performed 
using  room-temperature,  static  Ar  gas  of  known  density  and  polariza¬ 
bility,  a  method  of  calibration  or  normalization  was  employed  during 
the  course  of  the  measurements  which  obviated  corrections  caused  by 
position  changes  of  the  traversing  mechanism  of  the  source  incurred  by 
the  cooling  of  the  mechanism  from  room  temperature  to  its  much  lower 
steady- state  value.  This  procedure  is  mentioned  in  a  following  section. 

The  photomultiplier  tube  dark  count  rate,  determined  with  the  laser 
source  off,  and  the  background  count  rate,  determined  with  no  gas  within 
the  scattering  volume,  were  assumed  to  be  additive  corrections  and  were 
simply  subtracted  from  the  experimental,  photon  count  rate  values. 


3.0  EXPERIMENTAL  RESULTS  AND  ANALYSIS 


3.1  RESULTS 

Sonic  orifice  expansion  flow  fields  with  reservoir  pressures  of  200, 
250,  350,  500,  and  750  torr  were  produced,  and  Rayleigh  scattering 
measurements  were  performed  over  the  axial  distance  range  of  0.  4  ^ 
x/D  ^  25.  The  reservoir  temperature  employed  was  constant  for  each 
value  of  PQ  but  varied  from  275°K  at  250  torr  to  281°K  for  PQ  =  500  torr, 
and  the  data  are  reported  for  the  nominal  reservoir  temperature  of 
280°K.  Using  the  isentropic  flow  predictions  of  Sherman  and  Ashkenas 
(Ref.  12),  Fig.  3  shows  the  expansion  isentropes  plotted  in  the  P-T 
plane,  and  also  shown  is  the  vapor  pressure  curve  of  bulk  Ar  as  given 
by  Ref.  13.  The  crossing  of  the  vapor  pressure  curve  by  the  various 
isentropes  approximately  defines  the  values  of  P,  T  and,  thereby,  x/D 
of  the  flow  field  for  which  saturation  occurs.  Continuation  of  the  ex¬ 
pansion  along  the  isentrope,  once  having  crossed  the  vapor  pressure 
curve,  constitutes  a  supersaturated  flow  field  for  which  condensation 
onset  may  occur. 
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Figure  3.  Condensation  diagram  of  argon  for  a  sonic  orifice  expansion. 

The  experimental  results  are  shown  in  Figs.  4  through  8,  and  Fig. 

9  shows  all  results  on  a  single  graph.  The  scattering  data  were  normal¬ 
ized  at  one  position  for  each  value  of  PQ.  The  axial  position  selected 
had  to  be  sufficiently  far  from  the  source  that  flow  development  was  ful¬ 
filled,  and  the  equations  of  Ref.  12  were  valid  and  yet  sufficiently  close 
to  the  source  that  saturation  had  not  yet  occurred.  Also  of  consideration 
was  the  quality  and  reproducibility  of  the  data  at  the  point  in  question; 
i.e.,  it  was  obviously  desirable  to  select  the  best  data  point  available 
for  the  reference  value.  The  positions  selected  for  this  normalization 
are  indicated  in  Figs.  4  through  8.  Figure  7  shows  the  results  of  two 
separate  axial  scans,  and  the  scatter  shown  is  indicative  of  the  repro¬ 
ducibility  of  the  data.  It  should  be  noted  that  the  predictions  of  Sherman 
and  Ashkenas  for  the  axial  dependence  of  (Ni/N0)°  are  restricted  to 
x/D2  1,  and  the  experimental  data  shown  for  x/D  <  1  are  denoted  ex¬ 
plicitly.  It  is  interesting  to  look  more  closely  at  the  values  of  (Ni/N0)° 
shown  in  Fig.  9  for  the  region  x/D  <  1  and  consider  the  extrapolated 
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Figure  8.  Rayleigh  scattering  results  for  PQ  =  750  torr. 


Figure  9.  Axial  variation  of  effective  argon  density 
(N/N„)eff . 
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value  of  the  number  density  at  x/D  =  0.  The  initial  reaction  of  this  ex¬ 
trapolated  value  is  that  its  value  exceeds  expectations.  However,  Ash- 
kenas  and  Sherman  present  Mach  number  results  of  the  transonic  region 
for  a  sonic  orifice  flow  of  N2  and  show  that  at  x/D  =  0,  the  Mach  num¬ 
ber  is  approximately  0.  70,  not  unity.  This  value  of  Mach  number  cor¬ 
responds  to  an  isentropic  expansion  area  ratio  of  1.  0944.  Therefore, 
assuming  an  identical  area  ratio  at  x/  D  =  0  for  an  Ar  flow  field,  one 
finds  the  Mach  number  for  Ar  at  x/ D  =  0  to  be  0.  69  as  compared  with 
0.  70  for  N2.  Consequently,  at  x/ D  =  0  for  Ar  the  calculated  density 
ratio  is  approximately  0.  68,  and  it  is  seen  that  the  present  experimen¬ 
tal  results  are  not  inconsistent  with  this  value.  Moreover,  for  y  =  5/3, 
(Nj/N0)°  equals  0.  487  at  a  Mach  number  of  unity.  It  is  seen  from  Fig. 

9  that  these  present  results  indicate  the  development  of  sonic  flow  to 
occur  at  approximately  x/D  =  0.  40.  This  result  is  pertinent  for  all  ki¬ 
netic  rate  models  and  calculations  for  which  the  spatial  region  x/D  <  1 
is  of  importance.  As  an  example,  in  Ref.  2  calculations  of  the  number 
of  three-body  collisions  which  occur  in  a  sonic  orifice  flow  field  were 
performed.  The  assumption  that  M  =  1  at  x/ D  =  0  was  made,  and  a  func¬ 
tional  form  was  assumed  to  connect  M  =  1  at  x/D  =  0  and  M  =  2.  81  at 
x/D  =  1.  The  present  results  give  indication  that  the  lower  end-point 
of  M  is  in  error  by  approximately  30  percent,  and  it  is  noted  that  this 
error  occurs  in  the  higher  density,  higher  collision  frequency  region 
of  the  expansion. 

For  the  spatial  region  1  1  x/Di  (x/ D)q  it  is  seen  that  the  agree¬ 
ment  between  experimental  and  predicted  results  is  acceptable.  The 
deviation' of  the  measured  values  from  the  predicted  density  variation 
is  accepted  as  a  manifestation  of  the  growth  of  molecular  clusters  and 
the  onset  of  condensation.  It  is  noted  that  for  the  lowest  pressure  ob¬ 
served,  PQ  =  200  torr,  condensation  resulted.  This  result  is  to  be  con¬ 
trasted  with  the  previous  measurements  (Ref.  6)  of  N2  for  which  the 
isentropic  behavior  was  verified  over  the  axial  region  1  ^  x/D  ^  25  for 
PQ  =  500  torr.  The  axial  location  of  condensation  onset  (x/D)q  is  deter¬ 
mined  from  Figs.  4  through  8,  and  the  results  are  plotted  in  Fig.  10  as 
a  function  of  PQ.  From  Fig.  10  it  is  seen  that  significant  supersatura¬ 
tion  is  supported  by  Ar  for  PQ  of  approximately  200  to  250  torr,  the 
flow  having  proceeded  almost  three  orifice  diameters  beyond  saturation 
before  onset  occurs.  However,  as  PQ  increases,  ( x/D)q  decreases, 
and  condensation  follows  closely  on  the  heels  of  saturation  for  PQ  =  75  0 
torr.  In  fact,  the  results  indicate  that  for  reservoir  pressures  exceed¬ 
ing  750  torr,  condensation  can  be  expected  to  occur  at  saturation,  i.e., 
the  expansion  can  support  no  supersaturation  whatsoever.  The  super¬ 
saturation  ratio  s  =  Pg/P q,  the  ratio  of  saturation  and  onset  pressures. 
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^x^onset 

Figure  10.  Pressure  variation  of  condensation  onset  location  of  argon. 


is  easily  obtained  using  the  isentropic  predictions  of  Ref.  12  and  the 
results  of  Fig.  10,  and  s  is  shown  as  a  function  of  PQ  in  Fig.  11.  The 
experimentally  observed  degrees  of  supercooling,  Tg  -  Tg,  is  shown 
versus  reservoir  pressures  in  Fig.  12.  The  locus  of  the  onset  of  con¬ 
densation  in  the  P-T  plane  is  shown  by  the  dashed  line  in  Fig.  3,  and  it 
was  found  that  for  the  pressure  range  investigated  this  locus  was  ap¬ 
proximately  defined  by  the  relation 

P/T  „  lO'5-4 


22 


AE  D  C-T  R-74-32 


Figure  13  shows  the  axial  variation  of  the  depolarization  ratio,  p, 
for  PQ  =  75  0  torr,  and  it  is  noted  that  p  is  defined  as  the  ratio  of 
perpendicular-to-parallel  polarization  intensity  components;  i.e., 

p  - 

It  is  seen  by  comparison  of  Figs.  13  and  8  that  the  axial  region  of  p 
measurements  includes  the  initial  isentropic  portion  of  the  expansion 
as  well  as  those  of  condensation  onset  and  cluster  growth. 


Axial  Distance,  x/D 

Figure  13.  Axial  variation  of  depolarization  ratio  of  argon. 


3.2  ANALYSIS 

The  results  of  Figs.  4  through  8  are  used  with  Eq.  (6)  to  determine 
the  axial  variation  of  the  ratio  R,  and  Fig.  14  shows  the  result  for  the 
five  reservoir  pressures  investigated.  It  should  be  noted  that  for  x/D  = 
24  the  scattering  signal  for  PQ  =  200  torr  is  approximately  100  percent 
larger  than  its  isentropic  value,  and  for  PQ  =  750  torr  a  two-order-of- 
magnitude  increase  in  Is/l0  i-s  observed.  Of  further  interest  is  the 
steepening  of  the  growth  curve  as  PQ  increases. 

It  is  appropriate  to  consider  in  detail  at  this  point  the  argon  results 
of  Golomb  et  al.  (Ref.  2),  and,  in  particular,  focus  our  attention  on  Fig. 

6  of  this  reference.  It  was  found  in  this  work  that  the  dimer  peak  results 
scaled  according  to  PqD.  Table  1  presents  the  PQ  values  of  the  present 
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Figure  14.  Axial  variation  of  ratio  R. 


Table  1.  Condensate  Mass  Fraction  Results 


P0,  torr 

200 

250 

350 

500 

750 

PQ  (Refs.  2  and  5),  torr 

576 

734 

1008 

1440 

2159 

uc/ u°  (Ref.  5) 

1.  001 

1.  007 

1.  020 

1.  033 

1.  077 

g  (Ref.  5) 

<0.  001 

0.  015 

0.  042 

0.  07 

0.  16 

(N/Z)*(Ref.  8) 

20 

30 

55 

100 

200 

g  (this  work) 

0.  046 

0.  13 

0.  23 

0.  30 

0.  37 

Xc  (this  work)  x  103 

2.  4 

4.  8 

5.  5 

4.  4 

2.  9 
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investigation,  the  corresponding  PQ  values  of  Ref.  2,  as  well  as  the 
velocity  augmentation  and  condensation  mass  fraction  of  Ar  observed 
by  Bailey  (Ref.  5)  using  the  same  apparatus  but  in  a  later  work.  For 
the  200-torr  Rayleigh  scattering  results  of  Figs.  4  and  14,  it  is  seen 
that  at  the  reservoir  pressure  of  576  torr,  corresponding  to  the  orifice 
diameter  0.  39  mm  of  Refs.  2  and  5,  the  dimer  concentration  is  pre¬ 
dominant,  and  the  mass  fraction  of  condensate  is  on  the  order  of  0.  1 
percent  or  less.  Assuming  the  200-torr  results  of  Fig.  14  to  be  the 
result  of  dimer  scattering  only,  Eq.  (8)  yields  a  dimer  mole  fraction 
of  approximately  0.  25  and  a  condensate  mass  fraction  of  approximately 
0.  40,  which  differsjby  approximately  three  orders  of  magnitude  from 
Bailey* s  results.  If  one  attempts  to  take  account  of  the  obviously  exist¬ 
ent  trimers  shown  in  Fig.  6  of  Ref.  2,  it  is  found  that  X2/X3  «  6.  7 
from  that  figure.  It  is  easily  seen  that  Eq.  (6)  can  be  written  as 

where,  in  view  of  the  small  value  of  g  reported  by  Ref.  2  for  this  res¬ 
ervoir  condition,  the  approximation 

Nj  «  N2  «  Nt 

is  expected  to  be  quite  accurate.  Assuming  a^l  a\  «  i  and  solving  the 
above  equation  for  X2,  using  Eq.  (9),  one  finds 

g  «  0.32 

In  view  of  the  discrepancy  of  orders  of  magnitude  for  this  reservoir 
condition,  the  inclusion  of  the  trimers  is  not  significant.  Note -that  the 
results  of  Refs.  2  and  4  were  acquired  at  x/D  >  500,  whereas  those  of 
Fig.  14  span  the  range  of  1  <  x/D  <  25.  It  must  be  concluded  that  the 
present  results  indicate  scattering  from  clusters  much  larger  in  size 
than  either  dimers  or  trimers.  The  discrepancy  between  this  conclu¬ 
sion  and  the  results  of  Fig.  6  of  Ref.  2  can  be  explained  either  by  in¬ 
voking  a  difficult-to- accept  evaporation  process  of  the  cluster  between 
x/D  =  25  and  500  or  that  the  mass  spectrometric  results  previously  dis¬ 
cussed  quite  possibly  are  the  manifestation  of  fragmentation  processes 
of  larger  clusters. 

The  results  of  Hagena  and  Obert  (Ref.  8)  are  now  considered,  and 
it  is  desired  to  obtain  the  predicted  mean  value  (N/ Z)  appropriate  to 
our  values  of  PQ  and  D.  Reference  8  states  that  (N/Z)*  for  Ar  scales 
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approximately  as  P0D®'^,  and  Fig.  15  shows  the  approximate  results 
of  Ref.  8  for  the  variation  of  (N/Z)*  with  PQ  for  two  values  of  D. 
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Figure  15.  Variation  of  (N/Z)*  with  P0. 

Figure  16  shows  the  scaled  results  of  Ref.  8  for  a  3.  2-mm  orifice.  It 
is  seen  that  Hagena  and  Obert  predict  (N/Z)*  =  20  at  PQ  =  200  torr, 

D  =  3.  2  mm,  and  these  values  are  listed  in  Table  1.  Acknowledging 
both  the  distribution  of  cluster  sizes  encountered  by  Hagena  and  Obert 
and  the  absence  of' any  better  assumption,  the  authors  assumed  a  Poisson 
distribution  in  cluster  sizes  with  mean  value  J  =  (N/Z)  '  It  should  be 
noted  that  the  value  of  (N/Z)*  is  not  the  mean  value  even  if  the  distri¬ 
bution  were  Poisson;  rather,  it  is  such  a  value  that  one-half  the  ion 
cluster  current  is  caused  by  cluster  sizes  N/ Z  >  (N/Z)*.  Therefore, 
in  principle  one  should  plot  the  cumulative  Poisson  distribution  functions 


27 


AEDC-TR-74-32 


for  various  values  of  J  and  graphically  determine  the  best  value  of  J 
consistent  with  (N/Z)*.  The  Poisson  distribution  function  of  J  -  20  is 
shown  in  Fig.  17.  (N/Z)*  and  J  differ  by  approximately  two  percent, 

and  this  difference  is  henceforth  neglected.  Also,  it  should  be  noted 
that  for  J  =  20,  p(l)  <  10-^  and  is  negligible,  as  previously  asserted. 
Using  Eqs.  (19)  and  (21)  with  the  ascertained  values  of  R  and  the  values 
of  (N/Z)*  =  J  from  Ref.  8,  the  mole  and  mass  fraction  of  condensate 
Xc  and  g,  respectively,  are  calculated,  and  comparisons  are  made 
with  the  results  of  Ref.  5.  These  results  are  shown  in  Table  1.  It  is 
seen  that  the  g  results  differ  from  those  of  Bailey  by  over  an  order  of 
magnitude  for  PQ  <  250  torr,  but  the  agreement  improves  to  within  a 
factor  of  2  to  3  for  the  highest  pressure.  This  higher  pressure 
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discrepancy  is  perhaps  quite  acceptable  considering  the  crudity  of  the 
model  and  approximations  required  for  interpretation  of  Rayleigh  scat¬ 
tering  from  clusters  of  unknown  polarizabilities. 


Figure  17.  Poisson  distribution  function  for  a  mean  value  of  20. 


An  alternate  approach  is  to  assume  the  validity  of  the  values  of  g 
of  Ref.  5  and  to  calculate  J  from  the  experimental  values  of  R.  Table 
2  shows  these  results  and  those  of  Ref.  8  for  comparison.  Once  again, 
very  large  values  of  cluster  size  are  required  for  explanation  of  the 
Rayleigh -results.  However,  the  variation  of  J  with  PQ  is  somewhat 
peculiar,  and  the  procedure  employed  for  Table  1  calculations  is  pre¬ 
ferred.  Figure  18  shows  the  results  for  the  variation  of  g  with  P0. 

The  discrepancy  in  mass  fractions  obtained  by  Ref.  5  and  those  pro¬ 
vided  by  Rayleigh  scattering  for  the  low  values  of  PQ  is  not  understood. 
Although  it  is  possible  that  the  scattering  data  are  anomalously  high  for 
PQ  =  200  torr,  if  (N/Z)^  =  J  is  on  the  order  of  20,  one  intuitively  would 
expect  values  of  g  to  be  larger  than  0.  1  percent,  but  Bailey  reports 
them  to  be  less.  It  is  noted  that  the  flow  speed  is  a  measure  of  the 


Table  2.  Cluster  Size  Results  Using  Mass  Fraction  of  Ref.  5. 


PQ)  torr 

200 

250 

350 

500 

750 

g  (Ref.  5) 

£0.  001 

0.  015 

0.  042 

0.  07 

0.  16 

J 

1000 

300 

390 

590 

610 
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P0.  torr 

Figure  18.  Variation  of  mass  fraction  g  with  reservoir  pressure  P0. 


mass  fraction  of  condensate  so  long  as  certain  conditions  are  fulfilled, 
and  these  conditions  include  negligible  contributions  to  the  enthalpy  from 
both  the  translational  and  the  internal  modes  of  energy  of  both  phases  of 
matter.  Whether  violation  of  these  conditions  is  the  cause  of  at  least 
part  of  the  discrepancy  is  only  speculation  at  this  time. 

Turning  attention  to  the  validity  of  the  approximations  made  con¬ 
cerning  the  polarizability  cq  of  the  i-mer,  viz.  a-J a\  =  i,  this  assump¬ 
tion  was  made  as  a  result  of  ignorance  of  the  behavior  of  a  of  a  cluster. 
Obviously  the  absence  of  interaction  of  atoms  within  the  cluster  is  as¬ 
sumed,  or,  in  other  words,  the  effect  of  the  polarization  of  the  surround¬ 
ing  medium  is  neglected  when  calculating  the  electric  field  at  an  atomic 


l 
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site.  For  the  case  of  a  dimer  the  estimate  is  reasonably  good.  Recent 
work  of  Levine  (Ref.  14)  has  provided  an  approximate  correction  to  the 
assumption  a 2/^1  ”  2  which  amounts  to  only  three  parts  per  hundred. 

The  experimental  results  for  the  depolarization  ratio  p  shown  in 
Fig.  13  are  seen  to  increase  from  a  value  of  approximately  3  x  10"^  at 
x/D  -  1,  reach  an  initial  maximum  value  on  the  order  of  2  x  10"  2  at 
approximately  x/ D  =  2,  and  finally  decrease  to  1  x  10'^  at  x/D  =  25. 

As  discussed  previously,  an  increase  in  p  from  its  negligible  value  for 
the  monomer  Ar  manifests  the  existence  of  clusters.  It  is,  of  course, 
of  interest  to  compare  the  experimental  results  with  theoretical  pre¬ 
dictions.  Levine  and  Birnbaum  (Ref.  15)  have  considered  the  problem 
of  the  calculation  for  Ar2  of  the  polarizability  asymmetry  factor  j3, 
which  is  defined  in  terms  of  the  longitudinal,  ap,  and  transverse, 
polarizabilities  of  the  body  axes  of  the  cylindrically  symmetrical  mole¬ 
cule  as 


Z3  =  aL  ~  a-p 

They  assumed  the  interatomic  potential  of  Ar2  to  be  given  by  the 
Lennard- Jones  12-6  potential,  and  they  found  an  expression  for  j3  as  a 
function  of  the  interatomic  separation.  Assuming  their  nondimension- 
alized  distance  of  separation  to  be  on  the  order  of  unity,  for  the  pur¬ 
pose  of  order  of  magnitude  calculation,  the  expression  of  Ref.  16  for 
j3  gives  j3  =  0.  2.  Using  the  expression  (Ref.  11) 

p9  =  3/32/  (45  a2  +  3/3 2) 


one  finds 


P2  -  2.6  x  10'4 

It  is  seen  that  the  experimental  value  of  p  at  x/D  =  1  is  2.  9  x  10“^  and 
increases  upon  obvious  further  cluster  growth  as  one  proceeds  down¬ 
stream.  The  rather  good  agreement  between  the  measured  and  pre¬ 
dicted  values  of  the  depolarization  ratio  near  (x/D)p  makes  it  appear 
quite  likely  that  dimers  are  observed  during  the  earliest  part  of  the 
cluster  growth  process  but  that  larger  clusters  rapidly  assume  the  pre¬ 
dominant  cluster  specie  role  and  do  so  by  x/D  =  2  at  the  latest.  Since 
p  is  independent  of  density  in  the  first  approximation,  the  observed 
variation  of  the  depolarization  ratio  with  axial  distance  is  a  clear  indi¬ 
cation  of  a  changing  cluster  symmetry  and  cluster  number.  Referring 
to  Fig.  13,  it  is  seen  that  if  the  dimer  region  is  restricted  to  x/D  <j  2, 
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as  evidenced  by  the  p  data,  the  dimer  mole  fraction  X2  must  be  approxi¬ 
mately  less  than  or  equal  to  0.  025.  Consequently,  the  experimental  re¬ 
sults  for  PQ  =  750  torr  appear  to  place  a  maximum  value  for  X2  of 
0.  025  on  the  proposed  theoretical  descriptions. 

Regarding  comparisons  of  the  present  results  with  other  investiga¬ 
tions,  particularly  those  of  Yealland  et  al.  (Ref.  10),  it  was  not  possible 
to  accurately  employ  the  functional  relations  given  by  Ref.  10  because 
the  PQ  and  D  values  employed  in  the  present  study  lie  far  outside  the 
range  of  validity  of  that  work. 

Finally,  concerning  the  electron  diffraction  results  of  Ref.  4,  the 
orifice  diameter  was  0.  2  mm  and  the  experimental  data  for  PQ  of  2.  5 
to  3.  0  atm  at  room  temperature  indicated  a  short-range  order  of  the 
scatterer,  implying  a  liquid  phase  or  dimers  for  Ar.  If  one  assumes 
PqD  scaling,  the  pressure  value  corresponding  to  a  D  =  3  mm  is  ap¬ 
proximately  500  torr,  at  which  PQ  the  Rayleigh  scattering  results  are 
consistent  with  cluster  sizes  on  the  order  of  100.  This  result  is  an  ex¬ 
ample  of  the  limitation  of  diffraction  data;  assuming  short-range  order¬ 
ing,  or  a  quasi-liquid  phase,  information  concerning  the  size  of  the 
cluster  is  absent.  It  should  be  noted  that  the  smaller  orifice  of  Audit 
required  higher  values  of  PQ  than  employed  in  the  present  work  using 
a  3- mm  orifice  to  obtain  equivalent  condensate  results.  However,  the 
expansions  from  the  higher  values  of  PG  of  Audit  resulted  in  the  flow 
saturating  in  the  liquid  phase  regime  rather  than  the  solid  phase,  as  in 
the  present  case.  The  exact  nature  of  the  variations  to  be  expected  as 
a  result  of  this  difference  in  source  conditions  is  not  understood  at  this 
time.  Finally,  for  utilization  of  the  diffraction  technique  for  quantita¬ 
tive  interpretation  of  cluster  densities,  more  work  appears  to  be  re¬ 
quired  to  enable  the  correlation  of  scattered  electron  intensity  and 
specific  cluster  concentrations. 


4.0  DISCUSSION  AND  CONCLUSIONS 


The  Rayleigh  scattering  measurements  reported  herein  have  dem¬ 
onstrated  the  utility  of  this  technique  for  locating  condensation  onset  in 
a  sonic  orifice  expansion  of  argon.  Additionally,  this  location  provides 
information  concerning  the  supersaturation  ratio  as  a  function  of  reser¬ 
voir  pressure.  Cluster  growth  is  also  observed,  but,  of  course,  the 
data  are  ambiguous  in  the  sense  that  the  cluster  number  density  and 
polarizability  are  inseparable  without  additional  information.  Three 
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types  of  additional  information  were  employed  and  compared  for  con¬ 
sistency  with  our  results:  (1)  mass  spectrometric  detection  of  individ¬ 
ual  clusters,  (2)  time- of- flight  measurements  of  the  monomer  flow 
speed,  and  (3)  a  charge-normalized  mean  cluster  size  using  a  retarding 
potential  technique.  This  comparison,  as  has  been  mentioned  preiously, 
required  the  use  of  the  so-called  scaling  law  relating  some  character¬ 
istic  of  the  condensation  process  to  the  reservoir  pressure  and  orifice 
source  diameter.  It  was  found  that  the  Rayleigh  scattering  data  were 
irreconcilible  with  the  assumption  of  quantitative  accuracy  of  mass 
spectrometric  cluster  signals  of  Ref.  2;  this  result  was  not  shocking 
since  the  difficulties  of  faithfully  sampling  such  a  flow  are  well  known. 
However,  it  was  found  that  the  assumption  of  validity  of  the  cluster 
sizes  reported  by  Hagena  and  Obert  (Ref.  8)  did  enable  the  calculation 
of  reasonable  values  of  condensate  mole  and  mass  fractions.  These 
values  of  mass  fraction  were  in  order-of-magnitude  agreement  as  ob¬ 
tained  using  flow  speed  measurements,  except  for  the  lowest  pressures 
investigated.  For  these  no  quantitative  explanation  is  offered.  Addi¬ 
tionally,  the  measurement  of  the  depolarization  ratio,  p,  of  the  argon 
clusters  through  the  cluster  growth  region  provides  a  rather  interesting 
demonstration  of  the  spatial  symmetry  changes  which  occur  during  the 
accretion  process.  Moreover,  the  order-of-magnitude  agreement  be¬ 
tween  the  measured  and  predicted  values  of  p  over  a  finite  axial  dis¬ 
tance  appears  to  delimit  the  existence  region  of  the  dimer  in  the  direc¬ 
tion  of  the  flow  and,  thereby,  imposes  an  additional  condition  which 
must  be  fulfilled  by  any  proposed  theoretical  description  of  the  conden¬ 
sation  process. 

For  the  purpose  of  extracting  as  much  information  as  possible 
from  the  Rayleigh  scattering  data,  the  present  authors  have  resorted 
to  two  assumptions  which  are  worthy  of  note,  not  because  of  a  desire 
to  champion  their  validity,  but,  rather,  because  they  do  affect  the  con¬ 
clusions.  First,  a  Poisson  distribution  function  was  assumed  to  char¬ 
acterize  the  cluster  size  distribution.  Although  such  an  assumed  dis¬ 
tribution  is  admittedly  crude,  it  is  evidently  superior  to  a  Dirac  delta 
function  description,  and  for  small  mean  cluster  sizes  the  Poisson  dis¬ 
tribution  might  be  expected  to  be  reasonable.  Second,  the  polarizabil¬ 
ity  of  an  i-mer  has  been  assumed  to  be  i  times  that  of  the  monomer. 

The  basis  for  this  assumption  is  the  very  weak  binding  of  the  clusters, 
a  feature  which  will  help  to  maintain  the  independence  of  the  atomic 
species  and  thereby  insure  additivity  of  the  polarizability.  This  approxi¬ 
mation  is  obviously  quite  good  for  van  der  Waal’s  dimers,  but  its  ac¬ 
curacy  for  increasingly  larger  clusters  is  uncertain.  However,  despite 
any  inaccuracy  inherent  in  our  treatment,  an  inescapable  conclusion  is 
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the  caution  which  must  be  exercised  in  attempting  to  fit  various  kinetic 
rate  models  to  mass  spectrometric  dimer  signals  of  condensing- expan¬ 
sion  flow  fields,  such  as  has  been  done  by  Milne  et  al.  (Ref.  3).  The 
same  conclusion  obviously  does  not  apply  to  the  molecular,  or  effusion, 
source  work  of  Leckenby  and  Robbins  (Ref.  7),  because,  quite  simply, 
the  dimers  were  produced  by  equilibrium  processes  in  the  reservoir, 
the  molecular  beam  was  formed  without  collisions,  and  the  possibility 
of  fragmentation  effects  in  the  mass  spectrometer  are  minimal. 

Further  work  is  in  progress  which  is  designed  to  explore  the  vari¬ 
ous  scaling  laws  relating  reservoir  parameters  and  involving  other  gas¬ 
eous  species.  For  molecular  species,  rotational  Raman  scattering 
measurements  have  been  performed  through  the  condensation  onset  re¬ 
gion,  and  it  is  hoped  that  such  information  will  provide  additional  in¬ 
sight  for  the  understanding  of  two-phase  flow  fields. 
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NOMENCLATURE 

V 

D 

Orifice  diameter 

g 

Condensate  mass  fraction 

Iq 

Incident  laser  beam  intensity 

Is 

Scattered  laser  beam  intensity 

Isir  Isi 

Parallel  and  perpendicular  polarization  components  of  Is 

Optical  calibration  constant  for  Eq,  (1) 

N 

Number  density 

i 

Ni 

Number  density  of  i-mer 

j 

I 

j 

N/Z 

Number-to-charge  ratio  for  cluster 

(N/Z)* 

Mean  number-to- charge  ratio  for  distribution  of  clusters 

P 

Pressure 

v  t 

p(i) 

Probability  of  observing  i-mer  among  all  clusters 

R 

Scattered  intensity  ratio  defined  by  Eq.  (6) 

s 

Supersaturation  ratio 

T 

T  emperature 

tr 

Rotational  temperature 

- 

u 

Flow  speed 

X 

Mole  fraction 

Xi 

Mole  fraction  of  i-mer 

X 

Axial  distance 

a 

Polarizability  of  molecule 

ai 

Polarizability  of  i-mer 

«L>  “t 

Longitudinal  and  transverse  components  of  a 

13 

Molecular  asymmetry  parameter 

Au 

Change  in  u  attributable  to  condensation 

7 

Specific  heat  ratio 

36 
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AEDC-TR -74-32 

A  Wavelength 

p  Depolarization  ratio 

SUBSCRIPTS 

c  Variable  evaluated  for  condensing  flow 

5  Saturation  value 

T  Total  value 

o  Reservoir  value 

6  Condensation  onset  value 

SUPERSCRIPT 

o  Isentropic  value 

CGS  units  used  unless  otherwise  stated. 
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